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SUMMARY 


A literature survey of helium-II liquid-vapor phase separation is presented. 
Currently, two types of He-II phase separators are being investigated: porous, 

sintered metal plugs and the active phase-separator. The permeability Kp shows 
consistency in porous-plug geometric characterization. Passive phase-separation 
using porous plugs has been demonstrated for heat fluxes ranging from 1.0 mW/cm 2 to 
100 mW/cm 2 , with corresponding mass fluxes of from 0.1 mg/sec •cm 2 to 10 mg/sec*cm 2 , 
respectively. These results cover a Kp range of 10 cm to 10 cm . Both the 
heat and mass fluxes increase with K p . Downstream pressure regulation to adjust 
for varying heat loads and bath temperatures is possible. For large dynamic heat 
loads, the active phase-separator has shown a maximum heat-rejection rate of up to 
2 W and bath temperature stability of 0.1 mK. This report recommends that porous- 
plug phase-separation performance be investigated for application to SIRTF and, in 
particular, that plugs of 10 -9 cm 2 < Kp < 10~ 8 cm 2 in conjunction with downstream 
pressure regulation be studied. 


1 . INTRODUCTION 


There will be a number of difficult problems associated with the performance of 
future infrared space telescopes in the space environment. One crucial requirement 
is that of maintaining the IR detectors at temperatures less than 2 K while reject- 
ing sizable heat loads, and in so doing, maximizing lifetimes. Current cyrogenic 
technology utilizes the latent heat of liquid helium II for cooling. In such sys- 
tems, a porous plug may be used to define the position of the liquid-vapor interface 
in near zero-g fields. The feasibility of these approaches has been proved by the 
successful operation of the Infrared Astronomy Satellite (IRAS) (ref. 1), where heat 
loads and temperatures remained relatively constant. Similar programs have since 
followed the same path, for example, the Cosmic Background Explorer (COBE) and the 
Infrared Telescope (IRT) . The operation of future space missions, however, will 
impose more demanding requirements. In particular, the Space Infrared Telescope 
Facility (SIRTF) may require that variable heat loads ranging over a tenfold span be 
rejected, while IR detectors be maintained at 1.8 K and have operational lives of 
more than 2 yr. Therefore, there is a need to determine the limits of He-II phase 
separation as applied to SIRTF and other appropriate projects. 

Four U.S. groups have done recent investigations on porous-plug phase separa- 
tion. In addition, a German group has studied both porous plugs and the active 
phase-separator (APS) . Explanation and prediction of the performance of He-II phase 
separation has led to diverse and nonconciliatory efforts between the various inves- 
tigating groups. In short, each group proposes different explanations of the avail- 
able experimental data. Summarizing these abundant data in a single, comprehensive 
report is no simple task, because the variations between experimental systems and 
procedures are so wide ranging. Still, an attempt is made here to present the avail- 
able knowledge as a survey-type report that will provide future investigators with 
a quick and easy-to-use reference to the current published data and to the viewpoints 



of the various groups. In addition, areas recommended for further study are 
included. These recommendations will be directed toward the application of porous- 
plug phase separation to SIRTF . 

The report is divided into four sections. Section 1 reviews He-II hydrodynamics 
based on the two-fluid model and also presents theories relating the two-fluid model 
to phase separation. Sections 2 and 3 describe porous-plug and APS experimental 
systems and performance, respectively, and section 4 presents conclusions and 
recommendations . 


2 . THEORY 
Two -Fluid Model 

Liquid helium exhibits peculiar properties under certain conditions. It is 
these properties, among others, that qualify it for use as a cooling agent for IR 
detectors at temperatures below 2 K in a zero-g environment. 

Figure 1 shows the phase diagram of helium. There are two liquid phases of 
helium: He 1 and He II. He I obeys the laws of classical Newtonian mechanics; He II 

has the properties of superfluidity. The phase transition between He I and He II is 
referred to as a lambda transition and is characterized by the increase in the heat 
capacity to an infinite value at the temperature (T^) and pressure (P^) of the phase 
change, that is, the lambda point (inset of fig. 1). 

A peculiar property of He II is that a temperature gradient produces a pressure 
gradient. Figure 2 shows this experimentally. The name given to this phenomenon is 
the fountain effect, or the thermomechanical effect. It is exactly this phenomenon, 
together with an appropriate porous-plug geometry, that is used to contain the 
liquid in a zero-g environment for liquid-vapor phase separation. 

The model used to describe the mechanism of transport in He II has been termed 
the two-fluid model; it was independently proposed by both Tisza (ref. 2) and Landau 
(ref. 3). According to the two-fluid model. He II can be thought of as being com- 
posed of two separately interacting fluid components: a normal fluid component 

having finite viscosity and entropy and a superfluid component having zero viscosity 
and zero entropy. The combination of the two components thus gives rise to the bulk 
liquid properties, 

p = p + p (1) 

s n 

where p is the bulk-fluid density and P s /p 1 as T -* 0 K and P n /p 1 as 
T -> T^. Both the superfluid and normal fluid density ratios are functions of tem- 
perature. Figure 3 shows this experimentally measured relationship. The total 
laminar-flow description of He II, then, is described in terms of the independent 
equations of motion for the normal and superfluid (refs. 4-6): 

p dv /dt = -(p /p)VP - p SVT + p V 2 ^ (2) 

n n n s n n 

p d^ /dt = -(p /p)VP + p SVT (3) 

o o fc> 

From equation (3) for the superfluid component under steady-state conditions and 
laminar transport, the fountain effect is described according to 
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(4) 


VP T - pSVT 


The bulk mass-flux density j is just the vector sum of the two components 


pv = 


. *■>* 
p V + p V 

n n s s 


(5) 



"t 




( 6 ) 


Two major regions of heat transport are observed: laminar linear transport 

where q « AT, and nonlinear turbulent transport with q « (AT) 1 / 3 . To describe the 
laminar heat flow, the thermal energy content is carried by the normal fluid compo- 
nent (ref. 7): 


q = (pST)v n (7) 

Under conditions of zero net mass flow, where j = 0 and P n v n = -p s v s , the normal 
and superfluid components flow in opposite directions with respect to each other. 

That is, the normal fluid flows down the temperature gradient as a classical 
Newtonian fluid, while the superfluid flows up the temperature gradient producing 
the fountain pressure. Thus, transport is seen to be convective rather than conduc- 
tive. Experimental measurements of the apparent thermal conductivity (K app = -q/VT) 
shows a heat flux that can be up to 800 times greater than that of copper. 

For fully developed turbulent transport, commonly known as Gorter-Mellink (GM) 
turbulence (ref. 8), a mutual friction between the two components gives rise to the 
nonlinear heat-flux relation. 


q = P s ST(s|VT|/p n A GM ) 1 / 3 (8) 

where Aqm is the Gorter-Mellink parameter and is a function of pressure and tem- 
perature. Note that the GM region is independent of the flow duct diameter. 
Investigation of this transport has been extensive (refs. 8-12). 

Attempts to predict the critical transition between laminar and turbulent flow 
have incorporated the use of a dimensionless heat-flow Reynolds number (N^ e ) . 
Dimotakis (ref. 13) has proposed that for zero net mass flow, N^ e = (qlAg^/ST) where 
£ is the geometric characteristic length. The critical transition has been found 
to be NRe(c) = 1 f° r long, insulated cylindrical channels where £ = ltd. 


Phase Separation 

Numerous studies of He-II heat transport have been conducted under zero net 
mass-flow conditions (fig. 2). However, for phase separation, vapor is on one side 
of a thermo-osmotic filter (the phase separator) and liquid is on the other side. 

Figure 4 shows the general conditions required for satisfactory phase separation 
in a terrestrial environment. Under saturated vapor conditions, the downstream tem- 
perature T^ is less than the bath temperature T^ . Both T<j and T^ are less than 
the lambda temperature. This creates the necessary temperature gradient across the 
phase separator to produce the fountain pressure, AP^. It is the fountain pressure 
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that will contain the liquid helium within the Dewar. The liquid will stay confined 
only if AP t is greater than the sum of the applied and static pressure heads. 

Currently, two types of phase-separation devices are being explored: the active 

phase-separator (APS) and sintered porous metal plugs. The theoretical flow descrip- 
tion through both devices uses the two-fluid model equations that are referenced with 
respect to zero net mass-flow conditions. Phenomenological descriptions have also 
been published. In general, certain assumptions and corresponding equations are 
seen repeatedly throughout the literature; they are presented below. 

1. A heat balance around the phase separator, assuming saturated vapor condi- 
tions, requires that the heat rejected through the separator be proportional to the 
mass-flow rate. The proportionality constant is the latent heat of vaporization, X: 

Q = mX (9a) 

The heat-flux density is thus 

q = mX/A = jX = pvX (9b) 

2. The normal fluid assumes laminar Newtonian mechanics: 

VP = n V 2 v (10) 

n n ' 

the normal fluid viscosity being ri n . 

3. The steady-state equations of the two-fluid model are applicable to steady- 
state phase separation. 

4. (Specific to porous-plug phase separators.) The flow is assumed to be 
laminar Darcy flow with the permeability Kp a characteristic of the porous media 
flow geometry: 

V/A t = KpVP/q (11) 

where A-p is the superficial cross-sectional area of the plug, V is the volumetric 
flow rate, and q is the fluid dynamic viscosity. 

Active phase-separator - Theoretical transport in the APS has been studied by 
Schotte. In brief, the APS utilizes an annular gap to obtain phase separation. A 
more detailed description of APS design and operation is given in section 3. 

Schotte (refs. 14 and 15) and Schotte and Denner (ref. 16) have presented a 
theoretical model for He-II transport. In a moving coordinate system with velocity 
v, equation (7) becomes 


q = pST(v n - v) (12) 

Combining this with equation (9b) yields 

pvX = pST(v n - v) (13) 
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which may be solved for pv: 


pv = pSTv n /(ST + A) (14) 

The normal fluid assumes flow according to the Hagen-Poiseuille equation, 

v = (d 2 /12n )VP (15) 

n n 

Substitution of pv for the measurable superficial quantities m/A and of equa- 
tion (15) into equation (14) yields the total mass flow through the phase separator, 

m = Z(p/n ) [ST/(ST + A)]AP (16) 

where the geometry factor is Z = (2tr/12)Rd 3 /£ for an annular gap; R is the radius 
of the annular gap; and d is the channel width. Equation (16) reduces to 

m = Zp/n (ST/A) AP (17) 

at low T where A >> ST. Equation (17) can also be obtained for a stationary 
coordinate system (v = 0). 

Porous-plug phase separators - A number of groups have studied phase separation 
through porous plugs. The principal investigators of the various groups are 

1. M. DiPirro, currently associated with COBE 

2. T. H. K. Frederking et al., UCLA Cryogenic Engineering Laboratory 

3. J. Hendricks et al., associated with IRT 

4. D. Petrac and P. V. Mason, associated with IRAS and the Superfluid Helium 
Experiment (SFHE) 

In addition, studies have also been conducted by 

5. M. Murakami et al., Institute of Structural Engineering, University of 
Tsukaba, Japan 

6. H. D. Denner, G. Klipping et al., Fritz-Haber Institut der Max-Plank- 
Gesellschaf t , Berlin, and U. Schotte, Institut fur Theoretische Physik, Freie 
Universitat, Berlin, associated with the German Infrared Laboratory (GIRL) 

7. G. R. Karr and E. W. Urban, Marshall Space Flight Center (MSFC) , IRT project 

Following is a brief description of the theoretical and phenomenological equa- 
tions used by the various investigators in their studies of porous-plug phase sepa- 
rators. Where appropriate, equations are compared between groups. 

Hendricks and Karr (refs. 17 and 18) have proposed flow models based on a number 
of possible flow geometries. Poiseuille flow in straight channels of diameter d is 
assumed, with d determined in different ways. Schotte 's flow description (eq. (16)) 
is used with basic variations to the Z factor. Using 
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( 18 ) 


Z = (A / i) (d 2 /32) 

P P 

where Ap is the total pore area available for flow, three cases have been presented 

1. Ap = eAj, where e is the plug porosity 

2. A = N dji where N„ is the number of grains and d is the flow channel 

-..Poo o o 

diameter 

3. Ap = NpTTdp/4, where N p is the number of pores and dp is the pore diam- 
eter based on the filtration test 


Evaluations have been conducted to determine the validity of these models based 
on the prediction of the laminar-to-turbulent transition region (ref. 18). Results 
obtained indicate that the geometry of item (2) above best describes the geometry of 
porous plugs, but larger by a factor of 10 than that predicted under zero net mass- 
flow conditions (ref. 18). Opinions about the mechanism of the transition to turbu- 
lent flow seem to be of current interest. Many investigators (ref. 19; P. V. Mason 
1984: personal communication; J. B. Hendricks 1984: personal communication; and 

M. DiPirro 1984: personal communitation) believe that the liquid-vapor interface 

retreats within the pores of the plug during turbulent flow. The mechanism of trans- 
port is controversial, however. 

DiPirro et al. (ref. 19) and DiPirro (ref. 20) have considered the effects of 
the surface tension in determining the position of the interface with respect to the 
plug. A force balance over the interfacial surface under saturated vapor conditions 
results in 

pSAT - (dP/dT) AT - pgh = o(dA/dV) (19) 


where a(dA/dV) is the surface free-energy per unit area. For conditions in which 
the fountain effect becomes large, the interface Is thought to retreat within the 
plug until equation (19) is satisfied. Equation (19) predicts the point at which a 
steep rise in AT should be observed for increasing T^ . Data, as obtained by 
DiPirro et al. (ref. 19), show the transition to be within reason, as predicted by 
equation (19) . 


Karr and Urban (ref. 21) have proposed the description of He-II heat transfer in 
the form of Fourier's law. 
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-K VT 
app 


( 20 ) 


The apparent thermal conductivity K app is assumed to be dominated by the con- 
ductivity of the liquid. Taking the apparent thermal conductivity of the liquid 
(ref. 22) to be K app = p 2 S 2 Td 2 /12n n , and Aj scaled by the porosity, e, gives the 
mass flux through a porous plug, 

m/Ap = p(d 2 /12)e(AP/n n ) (ST/A) (21) 

For predicting the laminar-to-turbulent transition, the effort of Dimotakis 
(ref. 13) is used. Under conditions where the superfluid velocity is less than a 
critical value, v c , laminar flow is predicted by equation (22) : 

v /v < d 3 /£(TreAPSTA /12n A) (22) 

sc GM n 
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Yuan (ref. 23), Yuan and Frederking (ref. 24), and Lee et al. (ref. 25) have proposed 
normal fluid Darcy flow as a solution to the steady-state force balance between the 
normal and superfluid. Addition of equations (2) and (3), neglecting inertial 
effects, yields 

VP = n V 2 v (23) 

n n 


which is just a Newtonian description of the normal fluid, equation (10). A solution 
to equation (23) can be written in the form of normal-fluid Darcy transport subject 
to Newtonian boundary conditions: 


v 

n 


= K VP^/n 
pn T n 


(24) 


It is noted that in equation (24), VP = VP T (eq. (4)) and K pn is the normal-fluid 
permeability analogous to Kp of equation (11) . 

The nondimensional form of equation (24) is written as 


with 


\ " n vt 


N = qK 1 / 2 / n ST 
q pn n 

» VT " P ! S|vT|K=/ 2 /,4, 


(25) 


(26) 


where N p is a normalized heat-flow Reynolds number, and Nyp is the thermal driving 
force (refs. 24 and 26). Equation (23) applies to zero net mass-flow conditions. 
Geometric characterization of a porous plug at liquid helium temperatures has been 
found to be useful, using K pn as the characteristic dimension. For phase separa- 
tion with superfluid/normal counterflow, equations (6), (7), and (9b) combine to give 

|vj = (P n /P - ST/A)q/(p s ST) (27) 

At low temperature, where p n /p << ST/A, an asymptotic approach to zero net mass 
flow is predicted (refs. 25 and 27) and equation (25) applies. Figure 5 shows this 
approach to j = 0 for various porous media tested under phase-separation conditions 
(T. H. K. Frederking 1984: personal communication). 


For fully developed Gorter-Mellink convection in a long, insulated channel where 
the channel diameter d is the characteristic length (ref. 30), 


N q P/P S 


K GM^ N VT P s /|,p n^ ^ 


(28) 


Note that equation (27) is independent of d. Evaluation of data for zero net mass 
flow with NyT p s^ p n ranging over 10 orders of magnitude shows Kq^ = 11.3 ±1.4 
(K(^ is the universal Gorter-Mellink constant) . For smaller ducts with diameters 
of the order of 0.1 to 10 ym, as encountered with porous plugs, K pn shows possible 
duct diameter dependence (ref. 27). These smaller ducts correspond to permeabilities 
of the order of 10~ 9 cm 2 to 10 8 cm 2 . Equation (27) has been applied with relative 
success to quasi-steady phase separation (refs. 27 and 31). 
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Petrac and Mason (ref. 32) have given a phenomenological description of phase 
separation based on experimental data. The mass flow through the plug is found to 
be proportional to the applied pressure gradient under laminar conditions, the pro- 
portionality constant being related to the plug geometry, 

m = (FA)AP (29) 

where F is the plug geometry factor involving the "flow-channel size," fluid vis- 
cosity, and superficial plug thickness. Equation (29) is related to K pn of the 
thermo -osmotic Darcy equation by 

K = (F/pST)£n(AP/AP T ) (30) 

pn i 

The F-factors of a variety of different porous media tested by Petrac and Mason are 
included in table 1. 


3 . POROUS-PLUG INVESTIGATIONS 
Test Systems 

The general system employed for studying He-II phase separation using porous 
plugs is depicted in figure 6. The porous plug is bonded at the end of a thermally 
insulated vent line which, in turn, is dipped into a He-II bath. The He-II bath 
represents the Dewar system. System parameters that are measured are the mass-flow 
rate through the plug, m; the downstream and bath temperatures, T^ and T^, respec- 
tively; the downstream and upstream bath vapor pressures, P^ and P^, respectively; 
the liquid level above the plug, h; and the applied pressure difference across the 
plug, AP = P^ - P^ (downstream and upstream are referenced with respect to the finite 
mass flow of helium through the plug) . Operational conditions for phase separation 
usually require T^ < T^ < T^ and AP + pgh < APrp. However, under certain condi- 
tions, T^ may be greater than T^, provided that T^ < T^ (ref. 32). In this case, 
the lambda transition occurs within the plug. A heater in the bath may be used to 
simulate various heat loads on the Dewar system. Downstream vent-line impedance may 
or may not be adjustable. 

Preferred plug materials are low-conductivity metals and ceramics. The plug 
may be inserted in the vent line using an indium 0-ring seal and pressure coupling 
(refs. 17 and 19) or it can be bonded to the vent line by either epoxy adhesives 
or welding (ref. 27; J. B. Hendricks 1984: personal communication; M. DiPirro 1984: 

personal communication) . Temperatures are measured by germanium or carbon resistance 
thermometers. Downstream and upstream pump filters and liquid N 2 cold traps are used 
to minimize plug contamination and to ensure accurate mass-flow readings. Vent-line 
radiation shields should be incorporated where appropriate. 

When operating, heat absorbed by the bath will be rejected through the plug via 
the latent heat of the finite liquid mass flow. Ideally, the liquid-vapor interface 
occurs at the downstream surface of the plug, though it has been hypothesized that 
under certain conditions, the interface may retreat within the plug pores (ref. 19; 

J. B. Hendricks 1984: personal communication; M. DiPirro 1984; personal communica- 

tion) . Conditions may also result in liquid downstream of the plug (breakthrough 
conditions). Under ideal operation, laminar transport is realized. Consequently, 
there is minimal resistance to the heat flow and large m's are obtainable at 
smaller AP's. In this laminar region, dm/dAP may be of the order of 
10 mg/ (sec *Torr) (ref. 17), but this is dependent on K p . For turbulent flow. 
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dm/dAP may be reduced by two orders of magnitude, hence operation in the turbulent 
region places severe constraints on the maximum amount of heat that can be rejected 
by pure phase separation for a given plug area. Between laminar and turbulent flow 
is the transition region where dm/dAP decreases as AP increases; it can be char- 
acterized by a critical m c and critical AP C . 

a • 

When designing a system to test the maximum Q capabilities of a plug, Q ma x, 
pressure-drop considerations will be important. The maximum AP across the plug is 
determined by the difference between the bath vapor pressure and the total downstream 
vent-line pressure drop . This, in turn, dictates the limits on 1 %^ and Qmax* T° 
obtain higher m max and Q m ax> either a greater AP must be achieved, or the plug 
resistance must be decreased. The first alternative can be obtained by reducing the 
pressure drop downstream in the vent line. The second alternative is possible by 
altering the plug geometry: pore size, area, and thickness. It is noted that the 

maximum AP's available will be characteristic of the pumping system employed 
(refs . 27 and 33) . 

Experimental parameters commonly recorded are m, AP, T,j, T b , P(j, and P b . 

Since there is only one degree of freedom, the bath temperature T b is often fixed, 
and data are taken under steady-state conditions (refs. 17 and 34; D. Petrac 1979: 
unpublished report). A typical set of data is exhibited in figure 7(a) (D. Petrac: 
unpublished report). To obtain a different T b , a combination of heater input to the 
bath and of outer bath vapor pressure or downstream vapor pressure or both must be 
appropriately adjusted. In this manner, steady-state bath isotherms are obtained. 

Some investigators have studied phase separation under quasi-steady-state condi- 
tions (refs. 21, 29, 31, and 33). A representative set of data is given in fig- 
ure 7(b) (ref. 29). These results can be compared to constant bath temperature 
results only if dT b /dt is kept small or if transient effects are accounted for. 

In the latter case, the mass-flow rate is then m = pVC dT/dt + Q^ n , where C is 
the system heat capacity per unit mass, t is the time, and Q in is the total heat 
input to the Dewar (applied plus parasitic) . Dynamic response of a porous-plug 
system would be important, however, where transient Dewar heat loads or stored heat 
that needs to be rejected are of interest. Such would be the case when the vent line 
is opened after a launch to cool the system down to operating temperatures or when 
only a small amount of helium remains and the total system heat capacity cannot ade- 
quately maintain temperature stability. 


Flow Control 

Methods of flow control are of considerable interest. Many schemes of flow 
modulation have been proposed and investigated. In particular, flow control by use 
of a downstream pressure regulating device shows promise. This scheme has the advan- 
tage of simplicity of operation and of minimizing the number of cryogenic mechanical 
components . 

The use of a downstream heating element to vary the temperature difference 
across the plug has also been studied (refs. 21 and 34). In this system, the control 
of the fountain pressure is used to regulate the mass-flow rate. Figure 8 shows the 
change in m as a function of _ AT across a 10-ym ceramic porous plug with 
T b = 1.90 K and Am = m - m (Qext = 0), Qext d s the applied heat flow downstream 
of the plug (ref. 34). A linear relation between m and AT is observed, as pre- 
dicted by the two-fluid model (X is not a strong function of temperature in this 
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range of T) with Am inversely related to AT. The slope (dAm/dAT) = 4.2 mg/(sec’K) 
with AT ranging from 0 to 1.5 K. 

Figure 9 shows # the transient response of the downstream pressure as a function 
of step inputs in Q ext ( re f- 21). The response has been written as a time- 
dependent exponential function, with 

P d (t) = P x [l - a 1 exp (-b x t) ] (31) 

for increased step inputs in Q ex t an< 3 with P^^ the final pressure as t a ± 

the proportionality constant, and b x the time-constant. For decreased step inputs 
i n Qext 5 

P d (t) = P 2 [l + a 2 exp(-b 2 t)] (32) 

with P 2 , a 2 , and b 2 the final downstream pressure, proportionality constant, and 
decay time-constant, respectively. Table 2 shows the measured values for a 5-um 
ceramic plug : (ref. 21). No values for a-,^ and a 2 are given. 

■ Heating of the Dewar bath is also possible, with the design of the system skewed 
toward the maximum anticipated heat loads. An important disadvantage, however, is 
the wasteful boiling off of liquid at low heat inputs. 

A rotatable shutter assembly has been designed and tested (refs. 26, 27, and 31). 
In this system, the effective area available for flow through a plug is varied using 
a rotatable shutter assembly located on the downstream side of the plug. The shutter 
is driven by a motor via a connecting shaft. This proof-of-principle experiment has 
shown effective flow modulation of up to 60% from a mean value of 16 mg/sec. 

Other possibilities include a parallel plug arrangement in which each plug has 
its own vent line and value. If initially a number of valves are closed, then m 
can be increased by opening a value to a particular plug. The flow area would then 
be distinctly defined and the limiting restrictions to the maximum available AP 
could be reduced for laminar flow operation. However, liquid breakthrough may occur 
in the closed vent lines if a sufficient AT is not present, and thermo-acoustic 
oscillations may increase heat leaks. Backstream pressure regulation of this parallel 
plug-valve arrangement is also possible. 


Experimental Results 

Though a number of groups have investigated porous-plug phase separation of 
He II and though a large amount of information exists, there seems to be a wide dis- 
crepancy in the explanation of the experimental results. This section will not 
necessarily attempt to explain the abundance of data with respect to the various 
proposed theories. Instead, the available experimental data will be presented under 
a common frame of reference, when possible, in an effort to provide a systematic 
presentation of porous-plug phase-separation performance. Table 1 lists the various 
plugs and corresponding references to the literature. 

Plug characterization - It is common to characterize plugs according to the manu- 
facturer's stated nominal pore size or nominal filtration size; however, this type of 
characterization shows wide variation between manufacturers. Nonstandardized methods 
and procedures in determining these "nominal sizes" account for this confusion. In 
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particular, sizes may be determined by methods of particle filtration or retention, 
or by bubble pressure testing (refs. 36 and 37). Despite these inconsistencies, 
pore size is often used to distinguish plugs. This report will not characterize 
plugs by pore size, however. Instead, the permeability Kp, as defined by Darcy's 
law (eq. (11)) has shown a degree of consistency and thus will be used to charac- 
terize plugs. Still, pore size will be retained and referred to for reader 
convenience . 

Figure 10 shows Kp as a function of manufacturer's stated nominal size at 
room temperature (ref. 27). Mott Metallurgical Company and Pacific Sintered Metals 
(PSM) both use the particle filtration method in characterizing their plugs. The 
particle-retention size measurement has been conducted by independent testing facili- 
ties. A good relationship is observed between Kp and the nominal size for plugs 
of a single manufacturer, but there is a wide discrepancy between plugs of different 
manufacturers. Porous plugs from Mott show Kp to be an order of magnitude larger 
than Kp of PSM plugs for a given nominal size greater than 20 pm. At nominal 
sizes of 1 pm, Kp of Mott plugs is about 3 times as large as Kp of PSM plugs. 

The collimated plug of Petrac and Mason (ref. 32), which is a collection of finely 
drawn stainless steel capillary tubes (750,000, each with a diameter of 4 pm) shows 
a close correlation to the porous plugs of PSM. 

Figure 11 shows the experimentally measured relation between bubble-pressure 
pore size and particle removal size. No well-defined relation between these two 
parameters is apparent, although a relationship between particle removal size and 
pore size can be deduced. 

It has been found that Kp decreases with decreasing temperature (ref. 23). At 
liquid-helium temperatures under phase separation, it may be advantageous to charac- 
terize the plug according to the asymptotic limit to zero net mass flow, using Kp n 
of the thermo-osmotic Darcy equation as described in equation (24). In addition, 

Kp n has been found to approach the room temperature value of Kp as the zero net 
mass-flow asymptote is reached. Figure 12 shows this relation under phase separation 
for plugs of reference 29 and T. H. K. Frederking (1984: personal communication). 

Phase separation- For a given plug with a characteristic Kp, there is a unique 
steady-state relation between m, AP, and that will trace out a surface contour 

much like a three-dimensional T, P, and V phase diagram. Increasing Kp will 
result in a higher surface, with the possibility of extended transport in the laminar 
regime, all else being constant. From this steady-state diagram, one can project a 
plane onto a two-dimensional graph with either m, AP, or T^ as a constant parameter. 

A common way of showing data is on a graph of m versus AP for various con- 
stant T^'s. Figures 7(a) and 13 show this interdependence for the results of two 
studies (ref. 17; and D. Petrac 1979: unpublished report). Apparent is a large ini- 

tial slope (dm/dAP), corresponding to ideal laminar flow. With an increase in AP, 
the slope (dm/dAP) is drastically reduced. The point at which this transition from 
laminar to turbulent transport takes place is typically designated as m c and P c . It 
may be noted, however, that m c and P c characterize a region of transition rather 
than a distinct point. In addition, m c and P c increase with increasing Tfo . Thus, 
two distinct flow regimes have been identified, with the transition from laminar to 
turbulent transport occurring within some intermediate region. 

Laminar flow: Referring to figures 7(a) and 13, it can be seen that there is a 

severe restriction when operating porous-plug phase separation in the laminar-flow 
regime., That is, the transition to turbulence is encountered at very small AP's, 
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in the range of 0.24 to 1.5 Torr, depending on and plug K p . Figure 13 shows 

data for a 0.5-pm plug (ref. 17). Figure 7(a) shows terrestrial data for a 0.5-pm 
(Kp = 10 -9 -cm 2 ) plug that is currently flying in IRAS (D. Petrac 1979: unpublished 

report) . 

Figure 14 shows the m versus AP relation for various Tg for larger pore 
size and permeability plugs. Larger K p and larger Aj plugs can lead to an 
increase in m at substantially reduced values of AP. However, the information 
that has been compiled is not considered adequate to support definite conclusions. 

In particular, the data may not necessarily have been obtained under steady-state 
conditions at constant T^, . Transient heating owing to the system heat capacity may 
account for higher than expected m's with respect to steady-state T^'s. Further, 
the question of whether phase separation still exists or whether liquid breakthrough 
occurs at reduced AT's or large AP's is unanswered. 

Figure 15 shows a graph in which the superficial plug geometry is eliminated as 
a factor in system performance. The coordinates are plotted as mass-flux density 
(m/A' P ) versus VP = AP/Jt. The effect of decreasing K p results in a decrease in 
m/A p for a given VP, as one would expect. This inverse relation between m/Ap and 
K p suggests a Darcy-type approach in describing the superficial fluid flow. From 
an analysis such as this, the plug geometric parameters may be separated from the 
inherent fluid properties. If this is the case, then specification of a plug for a 
particular application can be done with thought to breakthrough conditions. For 
example, if a large amount of heat must be rejected and if the specifications for 
AP or AT across the plug create a liquid breakthrough condition, then a large-area 
plug of small Kp may be used. This, of course, depends on the ability to predict 
breakthrough conditions with respect to plug geometric parameters. Currently, there 
is a need to investigate breakthrough conditions in more detail. 

Transition region: The mechanism of the transition from laminar to turbulent 

flow for phase separation is currently a matter of uncertainty. Some investigators 
suggest superfluid critical velocities initiating Gorter-Mellink turbulence. Others 
propose the retraction of the liquid-vapor interface within the plug that results in 
transport through a high resistance vapor phase. Still others suggest a combination 
of the two. The prediction of the transition region is still under investigation. 

The actual determination must be made for each individual plug under testing 
conditions . 

Figure 16(a) presents results of the transition point of a 0.5-um plug (ref. 17) 
as a function of T^ . The critical mass-flow rate m c is seen to increase from 
1.7 mg/sec to 5.7 mg/sec with increasing T^ from T^ = 1.65 K to 2.1 K. After 
Tjj = 2.1 K, m c decreases to 4.8 mg/sec at T^ = 2.15 K. The decrease may be due to 
the high superfluid velocity initiating turbulence at this lower value of m as 

T b - t a . 

Turbulent transport: It is possible to operate a porous plug in the turbulent 

regime. An immediate consequence and restriction, however, is the quick approach to 
the maximum available AP as a result of the small values of dm/dAP. Figures 16(b) 
and 16(c) show dm/dAP versus T^ for both the laminar and turbulent flow regimes, 
respectively, of a 0.5-pm plug (ref. 17). 

From figure 16(b), dm/dAP for laminar flow ranges from 9 mg/ (sec ’Torr) at 
Tg = 1.65 K to 28 mg/ (sec *Torr) at T^ = 2.1 K. This is about 2 to 4 times greater 
than dm/dAP for turbulent flow at the same bath temperatures. This stronger m 
dependence on AP and T^ for laminar flow can give wider ranges of control for a 
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given phase-separation system. But it may be more difficult to stabilize and m 
in the laminar-flow region as opposed to the turbulent -flow region because of the 
greater dependence of m on AP. For plugs with larger K p , it is possible that 
laminar transport may be extended to higher m c 's and AP c 's with respect to plugs 
with smaller Kp . However, there is insufficient information about this effect. 

4. THE ACTIVE PHASE-SEPARATOR 


System Description 

The active phase-separator (APS) (refs. 38-40) is a dynamic device used to regu- 
late mass flow. Regulation of m compensates for varying heat loads and changing 
bath temperatures . 

Figure 17 shows a schematic of the APS. The APS, which is composed of a sliding 
pin situated within a bushing, is located at the exit vent line from the He-II Dewar. 

The defining geometry for fluid transport is thus an annular gap geometry of variable 

length. Tests have been conducted (ref. 38) using a parallel gap of 5 ym, a conical 
section with a gap variability of 5 to 15 ym, and an open circular groove. The diam- 
eter of the annulus is 20 mm, and the gap length £ varies from 0 to 15 mm. Under 

laminar flow conditions (with m < 7 mg/sec) , section AB of the pin is in the bushing 
Phase separation thus occurs within section AB. For increased mass-flow require- 
ments, the pin is extended out of the bushing to where point T (top of bushing) is 
within the limits of section BC of the moving pin. Here, turbulent flow is realized. 
For still further increases in m, section CD moves up until point T is within the 
section limits (CD) . 

Under conditions of large m, when section BC to CD is used, a liquid film of 
helium flows above the bushing top. The heat needed to vaporize this helium film is 
transferred from the outer bath to the inner helium film via conduction through the 
heat exchanger, H. In this way, large amounts of heat can be rejected (up to 1.5 W) 
without pure He-II convective heat transfer. Hence, pin design is important to the 
adaptability of the APS to specific applications. Section DE is required for pin 
guidance. The function of the ball closure, F, is to substantially reduce the mass 
flow when necessary, for example, during evacuation of a flooded heat exchanger after 
filling and topping off procedures (refs. 38 and 41). The ball closure is not leak- 
tight to superfluids, however. The entire pin assembly itself is driven by an elec- 
tromagnetic induction subassembly coupled to a bath temperature feedback control sys- 
tem. Experimental data have shown that m can be adjusted from near zero values to 
90 mg/sec at high bath temperatures, T^ = 2.1 K. Mass-flow rates up to 35 mg/sec 
have been noted at bath temperatures of 1.6 K. 


Experimental Results 

Figure 18 shows the dynamic response of the APS for step functions of the bath 
input. Stabilization of AT across the phase separator after steady state is of the 
order of 0.1 mK (ref. 40). For decreased step inputs in the bath heater, the direc- 
tion of the fountain effect is seen to be opposite to the desired direction, that is, 
a vector force acting in the downstream direction. This helps to reduce m and 
maintain a constant T^ . Evidently, at t = 6-8 min, a stabilized negative AT is 
observed, indicating a balance between the thermomechanical effect and the maximum 
closure position of the pin. The mass-flow rate m follows £ and AT, as might be 
expected, decreasing to approximately 10 mg/sec at maximum pin closure and a reversed 


13 



(negative) AT. No bath-temperature /time functions are given. It would be informa- 
tive to know the bath temperature stability and time-response curves. Figure 18(b) 
shows £ and AT for a step input in Q ex t at constant T b = 2.02 K and a constant 
AP = 13.9 mbar. The same performance is observed as shown in figure 18(a). In addi- 
tion, curves are indicated for gaseous He flow through the APS, while maintaining a 
stable bath temperature and stable AT. 

Figure 19(b) shows m versus the gap length £ for constant T^ and AP. A 
monotonic increase in m with decreasing £ is observed, with an observed propor- 
tionality to £ -1 , as predicted by equation (16). Figure 19(a) shows log m versus 
£. It can be seen that up to 100 mg/sec of helium can pass through the phase sepa- 
rator at 2.0 K. Also, m can be reduced to approximately 1.3 mg/sec at the same T^. 

Figures 20(a)-20(c) show m versus AP at different constant gap lengths, £ 
and T b . The solid lines indicate calculations based on equation (16). In the 
laminar-flow regime, data are predicted well. For higher turbulent flow, there is 
a large difference between what is predicted and what is observed, indicating the 
possibility of other flow phenomena besides GM convection. 

5 . CONCLUSIONS 

The conclusions of this survey can be summarized as follows: 

1. Porous plugs can satisfactorily function as passive helium II liquid-vapor 
phase separators. Heat fluxes of the order of 10 mW/cm 1 2 * (m/A T of the order of 

1 mg/sec) have been demonstrated for a 0.5-ym (K p = 10 -9 cm 2 ) stainless steel plug. 

2. Two regions of transport exist: a laminar region in which d(m/A T )/dAP is 

large and a turbulent region in which d (m/Ap) /dAP can be as much as 2 orders of 
magnitude less. 

3. The maximum heat flux is primarily limited by the maximum available pressure 
drop from the bath to a vacuum. It can be increased by enlarging the plug cross- 
sectional area A T or the permeability K p . The minimum q, however, will be 
increased, and subsequently breakthrough conditions must be considered with the 
latter . 

4. Downstream pressure regulation to adjust for varying heat loads and bath 
temperatures shows promise. The use of a downstream heater to control the thermo- 
mechanical force, hence, m, should not be dismissed as an alternative. 

5. The active phase-separator can reject a larger dynamic heat load at lower 

bath temperatures. Heat rejection rates of up to 700 mW at 1.6 K and 2 W at 2.1 K 

have been demonstrated. System pin and heat -exchanger design can be tailored to 
specific performance requirements. Moving components in a cryogenic environment must 
be considered. 

It is recommended that an in-house facility be designed and constructed to test 
the performance of porous-plug phase separation for application to SIRTF. Following 
are the specific areas for which more detailed knowledge is necessary. 

1. Testing of large plugs (10 -9 cm 2 < Kp < 10~ 8 cm 2 ; pore size range from 

2 to 10 ym) for the possibility of larger heat rejection rates and the study of the 

limitations imposed by breakthrough conditions. 
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2. The effects of superficial plug geometry, £ and A^. 

■ • 

3. Downstream pressure regulation for dynamic system control of T^, Q ext , and 
AP, and the design and implementation of a downstream feedback mechanism to maintain 
system stability. 

4. Overall correlation of all obtainable data on plug performance into a common 
frame of reference to predict plug performance. 
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TABLE 1.- POROUS-PLUG TEST DATA 


Plug 

number 

Plug and 
pore , a 
pm 

Nominal 

size.fr 

pm 

Kp 

T room 

cm 2 

Area, 

cm 2 

Thickness, 

cm 

e 

Material 

Refer- 

ence 

Bonding 

Experi- 

ment 0 

F-f actor , a 
g/(dyne *sec) 

Comments 

1 

M,f ,0.5 

3.92-b 

9.6X10 -10 

25 

0.61 


Stainless 

(e) 


PT 


(/) , (g) 

2 

M,f ,0.5 

3.92-b 

1.24xl0 -9 

25 

.61 


Stainless 

(e) 


PT 


(/),(£?) 

3 

M,f ,0.5 

3.92-b 

1.04x10”® 

25.3 

.61 


Stainless 

(e) 


PT 



4 

M,f ,0.5 

3.92-b 

1.03x10”® 

4.7 

.61 


Stainless 

(e) 


PT 


(g) , (h) 

5 

M,f ,0.5 

3.92-b 

1.07x10”® 

3.5/3.0 2 - 

.61 


Stainless 

(®) 

basd/ 

PS 


Final IRAS 













plugS 1 

6 

M,f ,0.5 



1.27 

0.64 


Stainless 

17 


PS 



7 

M, f ,2 



5.07 

.318 

0.31 

Stainless 

27 

Stycast 1266 

ZNMF 


M2511" 

8 

M, f ,2 



5.07 

.318 

.309 

Stainless 

27 


PS 


K2S2 k 

9 

M, f ,2 


5 .59 x 10” 9 

5.07 

.318 

.314 

Stainless 

27 

Stycast 1266 

PT 


ms3 k 

10 

M,f ,5 

25-30-b 

1.87x10”® 

25.3 

.64 


Stainless 

(e) 

RTV-180 

PT 


Original 













IRAS plug 

11 

M, f ,5 


3.1xl0” 8 

5.07 

0.318 

0.346 

Stainless 

27 

Stycast 1266 

PT 



12 

P,f ,2-5 


5.5xl0 _1 ° 

4.7 J 

.67 

.18 

Stainless 

28 


PS 



13 

P,f ,2-5 


5.59x10”® 

5.07 

.318 

.366 

Stainless 

27 

Stycast 1266 

PT 



14 

P,f ,5-15 


1.05xl0 -7 

5.07 

.636 

.36 

Bronze 

27 

Stycast 1266 




15 

P,f ,5-15 


8.5x10”® 




Bronze 

(m) 





16 

C,7-10 

7-10 

1.26x10”® 

20 

0.64 

0.61 

AI2O3 

33 


PT/PS 

1.7xl0” 7 


17 

B,4 


4 .9xl0” 9 (n) 

0.26 

.64 

3x10® (o) 

Stainless 

28 


PT/PS 

3.9xl0” 6 

(p) 

18 

H, 1 


lxlO” 11 

10.18 

.5 

0.4 

Al-Si 

29 


PS 



19 

S ,9-11 


3.4xl0- 9 

10.18 

.4 

.4 

Glass 

29 


PS 



20 

5 


5xl0” 10 

15.1 

.65 

.52 

A1 2 0 3 

21 

Epoxy 

PS 



21 

M 



12.97 

0.64/1.91 


Stainless 

19 


PS 


C0BE 

22 

Fab. 

2 

3.5xlO” 10 

1.96 

.15 


Stainless 

34 


PS 



23 

Fab . 

10 

1.2xl0“ 9 

3.80 

.7 


Ceramic A1 2 0 3 

34 


PS 



24 

Fab. 

10 

9.4xlO” 10 

3.80 

1.0 


Ceramic A1 2 0 3 

34 


PS 



25 

10 


3.77xl0” 9 

15.1 

.63 

0.384 

Nickel 

35 


PS 


(P) 

26 

10 


2.24x10”® 

15.1 

0.66 

0.677 

ai 2 o 3 

35 


PS 


(p) 

27 

0.5 


5.73X10” 11 

15.1 

.65 

.548 

ai 2 o 3 

35 


PS 


(P) 

28 

C 

<0 .5-b 


3.6 

.64 



32 



1.5x10”® 


29 

Py 36060 

2-2.5 


6 

.02 



32 



2.9xl0” 7 


30 

T,15fine 

<3 


1.8 

.64 



32 



1 . 1x10”® 



a M = Mott Metallurgical Co.; _f = filtration grade; PSM = Pacific Sintered Metals; C = Coors; = Brunswick; H = Haldenwanger (Berlin); 

J3 = Schott; Fab = fabricated by study investigator; = Pyrex; T = Tetraglass. "^Values measured by the investigator; suffix "b" indicates 
bubble pressure test. °PT = permeability test; PS = phase separation; ZNMF = zero net mass flow. “F-factor defined by equation (29). 
e D. Petrac 1979; unpublished report, /'phase separator for powers >0.3 W. ^Minimum steady-state cooling power 6-8 mW/cm 2 . "Phase separa- 
tion at 44 mW not possible. 2 '3.5 upstream, 3.0 downstream, owing to location of thermometer. PBASD = Ball Aerospace Division Specifica- 
tion BPS 11.08. kplug designation by investigator. ^From reference 15. m T. H. K. Frederking 1984: personal communication. n At low 
temperature. ''’Tubes per square centimeter. PCollimated capillary structure with 750,000 capillaries. PGraphs flow rate versus with 

T u and T,j at various heat inputs . 









TABLE 2.- TIME-CONSTANTS FOR A TRANSIENT 
HEATING EXPERIMENT - HEATERS LOCATED 
UPSTREAM OF AND DOWNSTREAM FROM A 5-ym 
POROUS PLUG [from ref. 21] 


Power, 

mW 

b i’ 
min 1 

Torr 

b 2i! 

nun 

^2 5 
Torr 

Liquid heater 

362 

0.295 

BW 


3.05 

285 

.287 

£ 

HI 

3.81 

Plug heater 

45.3 

0.33 

4.83 

0.479 

3.43 

101 

.412 

6.46 

.392 

3.33 

178 

.515 

8.97 

.476 

3.71 

276 

.603 

12.21 

.587 

3.78 

304 

.720 

13.04 

.637 

3.71 

381 

.769 

15.1 

.774 

3.14 

400.9 

.861 

16.0 

1.067 

3.43 

465 

.714 

17.1 

.509 

4.57 

Full power ~ 1000 

1.02 

32.3 

1.00 

4.57 
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Figure 1.- Phase diagram of helium. Inset: heat capacity as a function of 

temperature . 
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Figure 2.- Zero net mass flow: at steady state, the fountain pressure is checked by 

the pressure head, i.e., AP'p = pgh (n and s represent the flow directions of the 
normal fluid and super fluid components, respectively) . 



Figure 3.- Superfluid and normal fluid density ratio dependence on temperature. 




Figure 4.- General conditions for He-II phase separation: AP T > AP + pgh; 

T b < T d < T X . 
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O FREDERKING (24) 
5-15 /urn BRONZE 

Kp = 8.64 X 10 _8 cm 2 

□ PETRAC (28) 
COLLIMATED 
HOLE PLUG 

K p = 3.1 X 10“ 9 cm 2 

A KLIPPING (23) 

10 Mm GLASS 



N V j 


Figure 5.- Phase-separation data approaching conditions of zero net mass flow as Tj-, 
decreases; solid line represents the asymptotic limit to the thermo-osmotic Darcy 
equation, Nq = Nyp (T . H. K. Frederking 1984: personal communication). 
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TO DOWNSTREAM 
VACUUM PUMP 



Figure 6.- Schematic of a porous-plug testing apparatus for He-II vapor-liquid-phase 

separation. 
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0 “ 1 2 3 4 5 6 7 8 9 10 11 

AP, mmHg 


(a) m vs AP phase separation of in-flight IRAS porous plug: 0.5-yin Mott plug, 

Kp = 1 x io 9 cm 2 (D . Petrac 1979: unpublished report). 


Figure 7.- Measured parameters during phase-separation experiments. 







TIME, arbitrary units 


Figure 9.- Transient response of the downstream pressure when a step in power is 
applied to a heater located downstream from a porous plug (ref. 21). 


O MOTT, AIR 
A PSM, N 2 GAS 
□ PSM, WATER 
• PETRAC COLLIMATED 


PLUG (27) 



NOMINAL SIZE, jum 


Figure 10.- Room temperature permeability vs manufacturer's stated nominal size; 
Mott and PSM size measurements by filtration. 


PARTICLE REMOVAL SIZE RANGE, jum 


Figure 11.- Bubble pressure pore size vs particle removal size range for Mott and 

PSM porous plugs. 



Figure 12.- K p vs T b showing K approaching at the low 

P om _ Pn PT room 

temperature asymptote to zero net mass flow. 




AP, mmHg 


Figure 13.- Mass-flow rate vs pressure drop across a 0.5-ym stainless steel 

porous plug (ref. 17). 


O KLIPPING (23) 
V KARR (32) 

A PETRAC (30) 
O PETRAC (28) 
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Figure 14.- Mass-flow rate vs pressure drop for a variety of porous plugs: 
Kr. range = 10 -11 cm 2 to 10~ 9 cm 2 ; pore size range = 1 to 10 ym. 



O PETRAC4jum COLLIMATED (48) 

K p = 5 X 10 -9 cm 2 

O KLIPPING 10 fim GLASS (23) 

K p = 3.4 X 10" 9 cm 2 

A PETRAC 2-5 M m STAINLESS STEEL (28) 
K p = 5.5 X 10" 10 cm 2 


□ KLIPPING 1 jum AlSi (23) 
K p = 1 X 10 _11 cm 2 



Figure 15.- Mass-flux density vs pressure gradient for porous media with 

10 -11 cm 2 < Kp < 10 -9 cm 2 . 
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VAPOR 



Figure 17.- Schematic of the active phase separator. 
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< y- 



(a) = 1.602 K, AP = 3.7 mbar. 



(b) - 2.012 K; AP = 13.2 mbar with liquid in contact with APS; 

AP = 13.9 mbar with only cold helium gas flowing through APS. 

Figure 18.- Dynamic response of the active phase separator (ref. 21): AT and m 

as a function of l for step inputs in Q ext • 
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mg/s 




(a) Pin radius R = 10 mm, curves + and o; (b) Pin radius R = 10 mm. 

pin radius R = 5 mm, full lines. 

Figure 19.- Mass-flow rate vs £ for constant and AP of the APS (ref. 38). 
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